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ABSTRACT 

The  thermal  instability  of  incompressible,  rotating  fluid  saturated  in  a  porous  medium  is  critically  examined  both 
analytically  and  numerically  within  the  framework  of  linear  analysis..  The  analytical  discussion  provides  the  sufficient 
conditions  of  stability  and  instability  and  the  characterization  of  modes.  By  actually  calculating  the  root  of  eigenvalue 
equation  (  of  degree  3)  neutral  stability  curves  are  drawn.  The  numerical  results  show  the  effect  of  various  physical 
parameters  on  the  critical  wave  number  ac.  The  important  conclusion  is  that  medium  porosity  parameter  RD'  and  Prandtl 
number  Pr  show  a  dual  character  .  The  role  of  other  physical  parameters  such  as  rotation  parameter  T  ,  Richardson  number 
J  and  non-dimensional  parameter  R2  has  also  been  discussed  in  the  paper. 
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INTRODUCTION 

In  recent  years,  there  has  been  considerable  interest  among  geophysical  fluid  dynamicists  to  study  the  flows  and 
the  flow  instabilities  in  porous  medium  because  of  numerous  scientific  and  industrial  applications.  The  instability  of  fluid 
flows  in  a  porous  medium  under  varying  assumptions  has  been  well  summarized  by  Scheidegger  (1960)  and  Yih  (1961). 
While  investigating  the  flows  or  flow  instabilities  through  porous  media,  the  liquid  flow  has  been  assumed  to  be  governed 
by  Darcy's  Law  (1856)  by  most  of  the  research  workers,  which  neglects  the  inertial  forces  on  the  flow. 

Brinkman  (1947)  proposed  a  plausible  modification  to  Darcy's  law.  He  argued  that  the  equation  in  a  porous 
medium  must  reduce  to  the  ordinary  Navier-Stoke's  equation,  in  the  limit  of  zero  permeability.  This  led  to  the  celebrated 


Brinkman' s  model.  He  proposed  the  introduction  of  the  term 


f  \ 

E 


V  in  addition  to  —  ^2V  in  the  equations  of  fluid 
<P 


motion. 


Instability  of  compressible  or  incompressible  flow  has  been  studied  extensively  by  a  number  of  research  workers 
in  past  few  decades.  In  almost  all  such  investigations,  the  Boussinesq  approximation  is  used  to  simplify  the  equations  of 
motion.  In  an  important  paper,  Goel  et.al.  (1999)  have  examined  both  numerically  and  analytically  the  hydro  magnetic 
stability  of  an  unbounded  couple  stress  binary  fluid  mixture  having  vertical  temperature  and  concentration  gradients  with 
rotation. 

Bansal  and  Agrawal  (1999)  examined  the  shear  flow  instability  of  an  incompressible  visco-elastic  fluid  in  a 
porous  medium  in  the  presence  of  a  weak  magnetic  field. 
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The  important  results  obtained  include  the  extension  of  sufficient  condition  of  stabilities  of  Goel  et.al.  (1997)  to 
the  hydromagnetic  case,  characterization  of  modes  and  bounds  on  the  complex  wave  velocity  of  unstable  modes. 

In  this  paper,  An  attempt  has  been  made  to  study  the  thermal  instability  of  incompressible,  rotating  fluid  saturated 
in  a  porous  medium.  Following  Brinkman,  the  porosity  term  is  included  in  the  equation  of  motion. 

Important  analytical  results  obtained  in  this  paper  include  the  sufficient  conditions  for  stability  and  instability,  we 
have  actually  solved  the  eigen-value  equation  of  degree  three  obtaining  all  the  three  roots  which,  in  general,  are  complex. 
On  the  basis  of  this  data,  curves  have  been  plotted  to  provide  neutral  stability  curves  demonstrating  clearly  the  role  of 
various  physical  parameters. 

FORMULATION  OF  THE  PROBLEM 

An  incompressible,  heat  conducting,  viscous  fluid,  saturating  in  an  isotropic  porous  medium,  is  confined  in  the 
infinite  space  -°o  <  x,  y,  z  <  °°  having  z-axis  vertically  upward  so  that  the  gravity  acts  in  the  negative  z-direction  and  is 
subjected  to  uniform  rotation  £2  about  z-axis  chosen  in  the  direction  of  gravity.  Uniform  temperature  is  maintained 
along  z-axis. 

Following  Brinkman,  the  governing  equations  in  a  rectangular  frame  of  reference  (x,y,z)  rotating  with  fluid,  with 
porosity  corrections  are: 


dt  (f> 


=  -Vp  +  pg  +  fx 


v2- 


<f>  K 


v  +  —  (vxQ)  , 

0 


V.  v  =  0, 


—  +  -(v.V)T  =  kV2T 
dt  0 

and 

p  =  p0[l-a{T-TQ)\ 


(1) 

(2) 
(3) 


(4) 


Where  V  ,  p  ,  [1  ,  T,  K,  a,  and  kj  are  respectively  fuid  velocity,  density,  dynamic  viscosity  coefficient, 
temperature,  thermal  diffusivity,  thermal  expansion  coefficient,  medium  porosity,  and  medium  permeability,  g  =  (0,0, -g)  is 
gravitational  acceleration  and  the  suffix  zero  refers  to  the  values  at  the  reference  level  z=0. 

The  basic  state  for  the  problem  under  consideration  is  characterized  by 

v  =  (0,0,0), 
T  =  TQ-Pz,  (5) 
Q  =  {0,0,Q), 

Where  /?  may  be  either  positive  or  negative. 
The  condition  for  equilibrium  requires  that 
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dp 


+  pg=0  , 


dz 

p  =  p0[l+afiz] 

and  p  =  p0-p0g 
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z  +  -afiz2 
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(6) 
(7) 
(8) 


PERTURBATIONS  AND  NORMAL  MODE  ANALYSIS 

The  basic  state  characterized  by  eq.(8)  is  slightly  perturbed,  equations  in  perturbations  are  linearized  within  the 
framework  of  linear  analysis,  arbitrary  perturbation  f'(x,y,z,t)  is  decomposed  into  wave-like  components  as 


f(x,  y,  z,  t)  =  f  exp[o?  +  i(ax  +  by  +  czj] 


(9) 


where  a,  b  and  c  are  the  real  wave  numbers  and  a ,  a  time  constant,  is  complex,  in  general  and  after  eliminating 
various  physical  quantities,  the  following  eigen-value  equation  is  obtained  : 


(cr  +  Kl2) 


(10) 


Where  I2  =  a+b2+c2  and  m  =  a+b2. 
Introducing  the  transformations 

a  =  U°a   And (7,  m,c)  =  —(l*,  m",c* ), 
d  d 

We  get  the  following  eigen-value  equation  (after  dropping  the  asterisks),  in  non-dimensional  form  as 
a3  +Aa2  +Ba  +  C  =  0, 

WhereA  =  (l  +  2pr  )R212  +  2R~d, 


(11) 


(12) 


B  =  (2  +  pr)prR22l4+2R2(l  +  Pr)l2R-Dl  +  R-D2  +  ^--2' 
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It  is  important  to  note  that  the  values  of  J  can  be  both  positive  and  negative  depending  upon  whether  the 
temperature  decreases  or  increases  in  the  vertically  upward  direction. 


ANALYTICAL  DISCUSSIONS 

In  this  section,  we  shall  prove  some  important  results  with  the  help  of  equation 
Case-I:  C  <  0 

C  is  negative  under  the  condition 

Jm2(prR2l2  +  RD)  >   l2(  P2rRll6  +  2prR22R-D1l4  +  R2RDl2  +  R2Tc2 ) 


(12). 


(13) 


For  this  condition  to  be  true,  it  is  necessary  that  the  temperature  must  decrease  in  the  vertically  upward  direction. 
It  is  well  known  that  this  decrease  of  temperature  (J  >  0)  in  the  vertically  upward  direction  promotes  instability. 

Theoreml:-  System  is  unstable  under  the  condition  C  <  0. 

Proof:  -  If  C  <  0,  the  product  of  roots  of  the  equation  (12)  is  positive  so  that  at  least  one  positive  root  exists, 
implying,  thereby,  the  instability  of  the  system. 

Theorem  2:-  The  unstable  modes  which  exists  under  the  condition  C  <  0  are  non-oscillatory. 

Proof:  -  The  imaginary  part  of  equation  (12)  after  its  division  by  a  yields 


Oi 


207  + A- 


=  0 


Since  A  is  positive  and  in  view  of  theorem  -1,  modes  will  be  unstable  when  C  <  0.  Therefore  the  above  equation 
can  be  written  as 


2(jr  +  A  +  - 


C 


=  0 


(14) 


It  follows  that  Oi  must  necessarily  be  zero  so  that  the  unstable  modes  are  non-oscillatory. 

Theorem  -1  and  2  when  combined  together  yield  that  if  C  <  0,  then  the  modes  are  unstable  and  non-oscillatory. 

Case-II :  C  >  0 

Theorem  3:-  The  non-oscillatory  modes  are  stable,  if  B  >  0. 

Proof:  -  Equation  (12)  for  non-oscillatory  modes  (  at ■  =  0)  reduces  to 

a3r  +  Aa2r  +  Bar  +  C   =  0 

Since  A,  B  and  C  are  all  positive,  ar  should  be  negative,  implying,  thereby,  the  stability  of  non-oscillatory  modes, 
if  exist  in  the  system. 
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NUMERICAL  RESULTS  AND  DISCUSSIONS 

Equation  (12)  is  a  cubic  equation  in  cwith  real  coefficients  which  depend  upon  Rp,  Pr,  J,  R2,  T  and  the 
wave  numbers  a,  b  and  c  (through  I  and  rri).  We  have  calculated  all  the  three  roots  of  this  equation.  Our  aim  has  been  to 
examine  quantitatively  the  effect  of  various  parameters  such  as  the  Prandtl  number  Pn  Darcy-Reynolds  number  , 

Richardson  number  J,  rotation  parameter  T  and  the  parameter  R0  .  This  has  been  achieved  by  actually  calculating  the 
value  of  critical  wave  number  for  different  values  of  above  parameters. 

Figure  1:  Shows  that  dual  character  of  permeability  parameter  RD~' '.  For  J  =  5,  as  RD~'  increases  upto  2.75,  the 
range  of  stable  wave  numbers  decreases  and  the  range  of  stable  wave  numbers  increases  as  RD'  further  increases  beyond 
2.75.  This  figure  also  shows  a  destabilizing  character  of  J  as  for  a  fixed  value  of  RD',  the  critical  wave  number  ac  increases 
as  J  increases.  This  fact  is  also  evident  from  Figure  2,  3  and  5. 

Figure  2:  Shows  the  effect  of  Prandtl  number  Pr  on  critical  wave  number  ac  This  figure  shows  the  dual  character 
of  Prandtl  number  Pr.  For  J  =  5,  as  Pr  increases  from  0  to  3.25  approximately,  the  range  of  stable  wave  numbers  decreases, 
however,  the  range  of  stable  wave  numbers  increases  as  Pr  further  increases  beyond  3.50. 

Figure  3:  Shows  the  critical  wave  number  ac  for  different  values  of  rotation  parameter  T.  As  T  increases,  ac 
decreases  so  that  the  range  of  stable  wave  numbers  increases.  It  is  concluded  that  T  has  a  stabilizing  character  and  a  large 
value  of  Tis  required  to  stabilize  the  system  for  all  wave  numbers. 

Figure  4:  Shows  the  destabilizing  character  of  J.  As  J  increases,  ac  increases  which  decreases  the  range  of  stable 
wave  numbers. 

Stabilizing  character  of  non-dimensional  parameter  R2  is  exhibited  in  Fig. 5  .  As  R2  increases,  the  range  of  stable 
wave  numbers  increases. 

CONCLUSIONS 

The  analytical  discussion  provides  the  sufficient  conditions  of  stability  and  instability  and  the  characterization  of 
modes.  The  numerical  results  show  the  effect  of  various  physical  parameters  on  the  critical  wave  number  ac.  On  the  basis 
of  numerical  discussion  &  neutral  stability  curves  obtained  in  paper,  it  is  concluded  that  medium  porosity  parameter  RD~' 
and  Prandtl  number  Pr  show  a  dual  character.  The  rotation  parameter  T  and  non-dimensional  parameter  /?2,  which  depends 
upon  thermal  diffusivity/r,  have  stabilizing  character  and  Richardson  number  J  has  destabilizing  character. 
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APPENDICES 


Figure  1:  Critical  Wave  Number  Ac  Vs  RD 1  (R2  =  Pr=  0.5,  T  =  10) 
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